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A potential model for the interaction between a charmed baryon (Ac, Sc, and E*) 
and the nucleon (N) is constructed. The model contains a long-range meson (tt and cr) 
exchange part and a short-distance quark exchange part. The quark cluster model is used 
to evaluate the short-range repulsion and a monopole type form factor is introduced to 
the long-range potential to reflect the extended structure of hadrons. We determine the 
cutoff parameters in the form factors by fitting the NN scattering data with the same 
approach and we obtain four sets of parameters (a - d). The most attractive potential 
(d) leads to bound AcN states with = O'*" and 1+ once the channel couplings among 
Ac, Sc and E* are taken into account. One can also investigate many-body problems 
with the model. Here, we construct an effective AcN one-channel potential with the 
parameter set (d) and apply it to the 3-body AcNN system. The bound states with 
J = lf2 and 3/2 are predicted. 


1. Introduction 

Recent development of hadron spectroscopy revealed that there may exist various molecular 
bound states of hadrons (observed as hadron resonances) . In particular, the observation 
of the unexpected X, Y, and Z mesons and the follow by theoretical studies indicates that 
heavy quark molecules are more plausible. This can be understood from the balance between 
the kinetic term and the potential in the Hamiltonian: a heavier system has a smaller kinetic 
energy [1-3]. 

The above naive expectation motivates us to explore possible bound states composed of 
charmed (or bottomed) baryons yc(6) (^c(f))) ^c(6)) ■ • • ) the nucleon (N) or nucleus. This 
is, of course, a natural extension of the hypernucleus, which is a nuclear bound state with 
one or more strange baryons Y ( A, S, H ... ). Hypernuclear spectroscopy in the last decades 
played a key role in analysing structures of hypernuclei and extracting information on the 
YN and YY interactions. Because there is no two-body bound states in AN or T,N systems 
and it is difficult to perform direct scattering experiments for the hyperons, it is important 
to get information on their interactions from the three-body or heavier nucleus with strange 
baryon(s). 
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Channels 

1 2 3 4 5 6 7 

= 0+ 

= 1+ 

AcNCSq) EcNCSo) KNCDo) 

AcN{^Si) T,cN{^Si) T,*cN{^Si) AcN{^Di) 'ScN{^Di) 'S*N{^Di) T,*N{^Di) 


Table 1 The S-wave A^N states and the channels coupling to them [9] 


The idea of the charmed hypernucleus is, in fact, old [4-11]. It was pointed out [6] that 
the SU(4) symmetry, though it is badly broken, for the one boson exchange (OBE) models 
predicts a weaker attraction between Yc and N, because the K exchange is replaced by the 
D exchange, which is suppressed by the heavier mass of the D meson. Recent re-analysis 
with a model based on the heavy quark effective theories suggested the possibility of bound 
YcN states. However, the model has a difficulty in describing the short-range interaction 
and its prediction was very sensitive to the cutoff parameters in the one-boson exchange 
interaction [9]. 

In this paper, we construct a potential model consisting of the long-range one-pion and 
one-sigma exchange interaction and the quark exchange effect based on the quark cluster 
model for the baryon-baryon interaction [12, 13]. The latter provides a short range repulsion 
between Yc and N. Here Yc represents Ac, Sc and S*. The model also contains the off- 
diagonal components, which couple AcN, YcN and channels of the same total quantum 
numbers. By adjusting the meson exchange parameters to reproduce low-energy data of the 
NN system, we obtain four sets of the YcN potential. 

We solve the coupled channel Schrodinger equation using the Gaussian expansion method 
(GEM) [14] and obtain AcN bound states. The potential is also applied to the three-body 
AcNN system. Bound state solutions are found and their properties are analysed. 

The paper is organized as follows. In Sec. 2, we present our model of the YcN potential and 
explain how to determine the model parameters. In Sec. 3, the two-body YcN systems are 
analysed with the obtained potentials and the results are examined. In Sec. 4, we construct 
an effective one-channel AcN potential model and apply it to the three-body AcNN system. 
Gonclusions are given in Sec. 5. 


2. Model of YcN interactions 

We first consider the two-baryon systems, with isospin I = ^ and spin-parity J'^ = 0'*' or 1+. 
Table 1 shows possible channels contributed by AcN, YcN, and YcN for each J. The channels 
with the orbital angular momentum L = 0 and 2 will be coupled by tensor force, while those 
with the same L are coupled mainly by central force. We solve coupled channel Schrodinger 
equation in a hybrid potential model. The model includes a pion exchange potential and a 
scalar meson exchange potential for long-range interactions [9] and a short-range repulsive 
potential coming from the quark exchanges between the baryons [12]. 
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2.1. One boson exchange potential 

The meson exchange part contains four kinds of terms: spin-independent, spin-spin, spin- 
orbit, and tensor, 

~ (i, hi (^TT)-^Tr) ^) + -f^3(^7r)-^TT! j" 

{l)-jAYi{m^,A^,r) + {OLs)ij ^3("i^,A^,r) 

( 1 ) 

where i and j are the labels of the channels and CTr{i,j) and Ca{i,j) are the relevant coupling 
constants including the isospin factor. The spin dependent operators, Ogpin, 0>ten , and Ols, 
and their expectation values are given in Appendix A. The r-dependent functions Yi, 
and contain the cutoff parameters A^^ and Aq-. The explicit forms of the functions are 
given in Appendix B. We will determine the coupling strength and the cutoff parameters 
later. 


VaiiO) = Caii,j) 


nig 

Ibvr 


2.2. Short range repulsion from Quark Cluster Model 

In a previous approach [9], we considered exchanges of the vector mesons for the short 
range part of the l^N interaction. They, however, do not provide enough repulsion at short 
distances and result in very deep bound states with compact wave functions. As the wave 
function of the baryons overlap significantly at short distances, the quark exchange effect 
becomes important. We here employ the quark cluster model (QCM) for the short-range 
potential [12]. 

In QCM each baryon is a cluster determined made of three quarks. When the two clusters 
overlap at small r, the antisymmetrization among the (light) quarks induces a non-local 
interaction. When the two baryons overlap completely i.e., r = 0, all the six quarks occupy 
the lowest energy orbit with a single center. Such a state is approximately given by a product 
of the Gaussian wave functions with an appropriate symmetry according to the quantum 
numbers. The strength of the repulsive potential at r = 0 (Vq) is thus determined by the dif¬ 
ference between the energy of the single-centered six-quark states and that of the individual 
baryons. 

Go ~ {HH\H - (2) 

It turns out [12] that the strengths are sensitive to and determined dominantly by the 
color-magnetic interaction (CMI) given by 

VcM = -S • crj)(Aj • \j) (3) 

i<j 


where cr, and Aj are the spin and color operators of the i-th quark. The expectation values 
of the color-magnetic operator for a five light quark systems can be computed by using the 
formula. 


{Vcm) = S 


8N + -S{S + 1) + 2^2 [5C/(3)c] - 4C2 [SU{6)c 

O 


(4) 
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where N is total number of light quarks, S is the total spin, and C 2 [SU ( 5 )] is the quadratic 
Casimir operator, whose value is specified by the Young diagram [fi, ■ ■ ■ , fg], 


C 2 [ 517 ( 5 )] 


1 

2 


fiifi — 2z + 5 + 1) 


9 


(5) 


The overall strength /3 can be determined by the A(1232) - nucleon {N) mass splitting, 
which comes also from Vcm- From {Vcm)a ~ {^cm)n — ~16/3 = 293 MeV (exp.), we obtain 
/3 = 18.2 MeV. Table 2 shows the values of Vcm for various two-baryon (six-quark) states. 
These values are evaluated in the heavy quark limit, so that the heavy quark spin does not 
contribute. 


System 

Vo [MeV] 

(iViV)f=? 

450 

(iVV)f=i 

350 

(A,V)f= 0/2 

300 

{AcNf=l/, 

300 


/ ^cN 


\ 5=0 

/ 1 = 1/2 


100 0 \ / SciV \ ^ ^ 

0 oj’ 


166.7 -24.0 

-24.0 108.3 


Table 2 Expectation values in MeV of the color magnetic interaction (CMI) for the 
relevant channels in the heavy quark limit. For the correlated channels of ScV and 
the values are given in the matrix form. 


We assume, for simplicity, that the radial dependence of the QCM potential is given by a 
Gaussian, 

= ( 6 ) 

The range parameter b is supposed to coincide with the extension of the quark wave functions 
in the baryon. According to [13], the typical values for the NN interaction are about 0.54 ~ 
0.58 fm. For the YcN systems, we use two typical values, b = 0.5 and 0.6 fm, and compute 
the results. 

2.3. Coulomb potential 

We also include the Coulomb potential between the charged baryons. In this calculation, we 
consider only the two-baryon channels coupling to AcN{I = ^). For the AcU system (total 
charge Q = -|-1), there is no Coulomb effect. For the Q = 2 channels, A^p, and 
we assume that the mixing of F+p and is not modified by the Coulomb interaction 

and make an approximation for the Coulomb potential with the combination, 
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Then the ’’effective” Coulomb potentials for the YcN channels are given by 




AeAf 


(r) = 


ahc 


coulomb 

coulomb^ / Sr 


1 ahc 


K 


t.:n 


(r) = - 


1 ahc 


coulomb^' ^ 3 r ^ 

where a ~ 1/137 is the hne structure constant. 


( 8 ) 


2.4- Determining the potential parameters 

In the previous study using the OBE potential [9], we found that the results are very sensitive 
to the cutoff parameters. To remedy this problem, in the present approach, we adjust the 
parameters of the potential to reproduce the NN interaction data using the same model. In 
doing so, we fix the vr-baryon coupling constants and the short-range potential derived from 
the quark model. Then the cutoffs, and Ao-, and the sigma coupling constant, Cg-, are 
the parameters to be determined in the NN system. After that the results are generalized 
to the YcN system. In this method, we assume that the light mesons couple only to the light 
quarks and the cutoff parameters are common to NN and YcN systems. The details are 
given in Appendix B. 

In searching the three undetermined parameters, A,r) A^- and Co-, we use the following 
experimental values of the NN interaction, 

NN{^Sq) scattering length = —23.7fm, 

NNt^Si) binding energy(deuteron) = 2.22MeV. (9) 

To get appropriate values, we restrict the possible range of parameters to be 

= 500 ~ 900MeV, 

Ao = 900 ~ 1200MeV, 

C„{NN) = -64.0 ~ -225.0. (10) 

Unfortunately, we could not find a solution. Therefore we relax the condition that Ca is 
independent of and take Ca{J'^ = 0+) and Ca{J'^ = 1^) independently. Then we find 
solutions, which are given in Tables 3 and 4. In Table 3, we show parameters with a fixed 
A-n- = 750 MeV and varied A^- while Table 4 lists the results with varied A.^- for a fixed 
Act = 1000 MeV. These two Tables also present the differences between = O'*') and 

Ca{J^ = 1 +), 

Ac7 = |Cct(J^ = 0+)-C'ct(J^ = 1+)|. (11) 

From the results in the Tables, one finds that Ac is insensitive to the change of Aq-, but Ac 
increases rapidly for increasing A^^. The features for b = 0.5 fm and 0.6 fm are qualitatively 
similar. Since the QCM repulsion for h = 0.5 fm is weaker than that for b = 0.6 fm, the 
resulting is also smaller and Ac tends to be smaller for b = 0.5 fm. We suppose that the 
spin dependence of is small, so that a small Ac is favoured. 

The resulting NN potentials for = 0+ with the obtained parameters are plotted in 
Figs. 1 and 2. We find that some of the potentials (for small At^) are strongly attractive 
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6 

= 0.6[fm] 


b 

= 0.5[fm] 


Act [MeV] 

^(0+) 

c.(i+) 

Ac 

Cct(0+) 

a(l+) 

Ac 

900 

-213.16 

-161.29 

51.87 

-176.89 

-132.25 

44.64 

950 

-179.56 

-134.56 

45.0 

-148.84 

-110.25 

38.59 

1000 

-156.25 

-118.81 

37.44 

-129.96 

-96.04 

33.92 

1050 

-141.61 

-106.09 

35.52 

-116.64 

-86.49 

30.15 

1100 

-127.69 

-96.04 

31.65 

-106.09 

-79.21 

26.88 

1150 

-116.64 

-88.36 

28.28 

-98.01 

-72.25 

25.76 

1200 

-108.16 

-82.81 

25.35 

-90.25 

-67.24 

23.01 

Table 3 The NN 2-body parameters for A 

= 750 MeV 


6 

= 0.6[fm] 


b 

= 0.5[fm] 


A^[MeV] 

Cct(0+) 

c.(i+) 

Ac 

Cct(0+) 

a(l+) 

Ac 

500 

-148.84 

-139.24 

9.6 

-118.81 

-121.0 

2.19 

550 

-151.29 

-136.89 

14.4 

-123.21 

-116.64 

6.57 

600 

-153.76 

-134.56 

19.2 

-125.44 

-114.49 

10.95 

650 

-153.76 

-129.96 

23.8 

-127.69 

-110.25 

17.44 

700 

-156.25 

-125.44 

30.81 

-127.69 

-104.04 

23.65 

750 

-156.25 

-118.81 

37.44 

-129.96 

-96.04 

33.92 

800 

-156.25 

-110.25 

48.51 

-129.96 

-88.36 

41.6 

850 

-156.25 

-102.01 

54.24 

-132.25 

-79.21 

53.04 

900 

-156.25 

-92.16 

66.6 

-132.25 

-70.56 

61.69 


Table 4 The NN 2-body parameters for A^- = 1000 MeV 



Fig. 1 NN potential {J'^ = 0+) for 6=0.6fm Fig. 2 NN potential {J'^ = 0+) for 6=0.6fm 
and Act = 1000 MeV. and A.^- = 750 MeV. 


at short distances and may not be appropriate for the current purpose, although they all 
reproduce the ^Sq scattering length and the deuteron binding energy. 

Therefore we choose the parameters with A.^- = 750 MeV and Act = 1000 MeV and employ 
four sets of the most realistic potential parameters, given in Table 5. We call these models 
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“CTNN” potentials, as these parameter sets correspond to the NN experimental data. The 
YcN potentials derived for J'^ = 0"*“ are shown in Figs. 3 - 8. In Appendix C, we show the 
YcN potentials for J'" = and also the individual components of the CTNN-a potential. 



a 

b[fm] 

parameter a 

-67.58 

0.6 

parameter b 

-77.5 

0.6 

parameter c 

-60.76 

0.5 

parameter d 

-70.68 

0.5 


Table 5 The CTNN potential parameters 


3. YcN bound states 

In solving the Schrodinger equation for the two-baryon system, we employ a variational 
method with Gaussian trial functions, the Gaussian Expansion Method [14], where the radial 
wave functions are expanded by the basis states given by Gaussian functions with varied 
range parameters. 


f^max 

^lm{r) = Cnl^nlmir) 

n=l 

<Pnlmir) = (pnlir)Ylmir) 
(pniir) = 

1 1 

Tn ria” ^ 


( 12 ) 


This method has been applied to various bound state problems and was proved to give an 
accurate approximation to the eigenstate-energies and wave functions. 

The Y(.N binding energies are given in Tables 6 and 7. The A^n system = O'*' or 1+) 
without the Goulomb potential has a bound state for the YcN-CTNN parameter b, c, and d. 
The CTNN-b and c potentials allow very shallow bound states and the CTNN-d potential 
gives a larger binding energy. The result indicates that there is a strong QGM repulsion at 
short distances for the GTNN-a and b parameters. On the other hand, the AcP system with 
the Goulomb potential has a bound state only for the parameter GTNN-d. 

As a cross check, we also calculate the AcY scattering lengths for the GTNN potentials. 
These results are consistent with the binding energies (see Tables 6 and 7.) 

The above results show that the binding energies of the = 1"*“ states are larger than 
those of = 0'*'. This is consistent with the previous study [9]. 

Table 7 also shows the probabilities of the coupled channels. For all the parameter sets 
and the total angular momenta, the probability of AcN is more than 99 % . In the case of 
= 0+, the probabilities of (^Dq) are larger than that of ScN (^^o). On the other 
hand, in the case of = I"*", the respective total probabilities of T,cN and Ti^N are almost 
equal. Looking closely, T,cN{^Si) and T,*N{^Di) contribute largely to the results. These 
observations can be explained with the contributions of the strong tensor force between the 
A^N {L = 0) and / Y*N {L = 2) channels. 
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Fig. 3 YcN-CTNN potential for the AcN Fig. 4 YcN-CTNN potential for the ScA^ 
single channel. single channel. 



Fig. 5 FcA^-CTNN potential for the T,*N Fig. 6 YcN-CTNN potential for the AcN- 
single channel. A^cN channels. 



Fig. 7 YcN-CTNN potential for the AcN- Fig. 8 YcN-CTNN potential for the T,cN- 
T,*N channels. T,*N channels. 


4. AcNN systems 

4.I. Effects of channel couplings 

In the previous section, we find that the probabilities of the AcN component are almost 
100% and the mixings of T,cN or Ti*N are small. However, effects of the T,cN — Ti*N channel 
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= 0+ 

CTNN-a 

CTNN-b 

CTNN-c 

CTNN-d 

B.E. [MeV] 

- 

- 

1.72 xl0"3 

1.37 

(+ Coulomb) 




(0.56) 

scattering length [fm] 

-3.64 

-65.15 

130.93 

5.31 

probability(Ac V) [%] 

- 

- 

99.97 

99.29 

probability(Sc V) [%] 

- 

- 

7.0 xlO-3 

0.20 

probability(S*V)[%] 

- 

- 

2.1 xlO-2 

0.51 


Table 6 Binding energies of AcN (J^ = 0+) for the CTNN potentials. The probabilities 
of the TicN and Ti*N channels are also shown. 


= 1+ 

CTNN-a 

CTNN-b 

CTNN-c 

CTNN-d 

B.E. [MeV] 

- 

1.67 xlO"'^ 

1.91 xlO-2 

1.56 

(-1- Coulomb) 




(0.72) 

scattering length [fm] 

-4.11 

337.53 

39.27 

5.01 

probability(Ac V) [%] 

- 

99.99 

99.90 

99.23 

probability(Ec V) [%] 

- 

4.9 xl0"3 

4.9 xl0"2 

0.39 

(D-wave {^Di)) 

- 

4.5 xlO"^ 

4.6 xl0"2 

0.35 

probability (S* N ) [%] 

- 

4.6 xlO-3 

4.6 xlO-2 

0.38 

(D-wave (^Di)) 

- 

3.1 xl0"3 

3.2 xl0"2 

0.25 


Table 7 Binding energies of AcN (J^ = 1+) for the CTNN potentials. The probabilities 
of the coupled D-wave AcN and J^cN and 'S*N channels are also shown. 



AcN- 

T^cN - T.*cN 

Ac 

N 

AcN- 

-^cN 

AcN- 



0 + 

1 + 

0 + 

1 + 

0 + 

1 + 

0 + 

1 + 

CTNN-a 

-3.63 

-4.10 

-1.11 

-1.11 

-1.16 

-2.07 

-3.13 

-2.09 

CTNN-b 

-63.25 

398.67 

-2.62 

-2.62 

-2.78 

-6.74 

-20.84 

-7.00 

CTNN-c 

139.07 

39.96 

-3.01 

-3.01 

-3.19 

-8.61 

-48.56 

-9.00 

CTNN-d 

5.32 

5.02 

-28.59 

-28.59 

-44.65 

9.79 

6.01 

9.36 

B.E. 

1.37 

1.56 

- 

- 

- 

0.36 

1.04 

0.39 


Table 8 Scattering lengths in fm, and two-body binding energies in MeV for the cases 
that the coupled channels are reduced to AcN only, AcN-TicN or AcN-T,*^ 


coupling are important in binding the AcN system, because the single channel calculations 
of AcN do not show any binding solutions for all the CTNN-a ~ d potentials. 

In Table 8, we show the scattering lengths obtained for partially coupled systems, i.e., 
AcN, AcN — TicN, and AcN — It is found that only the CTNN-d potential allows 

bound states in the AcN — (O'*' and I'*') and the AcN — T,cN (1+) systems. In particular, 

effect of the (O'*") channel is large because no bound state is found without this channel. 
This indicates that the tensor force from the one-pion exchange, which induces the coupling 
between {^Dq) and AcN (^^o), is significant. Similarly, the contributions of the tensor 
force and T,cN — T,*N couplings are found to be important for the = 1'*' state. 
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Fig. 9 A cN effective potentials for J'^ = 0"*“ and 1"*“. 


4-2. Effective AcN potential 

In order to apply the obtained potentials to many-body systems with Ac, we construct 
effective single-channel AcN potentials, in which the effects of the T,cN — Ti^N couplings are 
absorbed in phenomenological parameters. We assume two-range Gaussian forms, 

+ Vf+e~i. 

(13) 

Here, bi and 62 are range parameters, and and are the strength parameters. They 
are determined from the CTNN potentials. For simplicity, we choose the same bi and 62 for 
J'^ = 0"*“ and !“*■, and rewrite Eq. (13) to be 

2 2 

= [K' + • crVf] + [K' + g-% ^ ( 14 ) 


(15) 


Vf =l{V^+ + 3Vl+), 

V4 

Now, we use only the CTNN-d YcN potential to construct the effective one. By reproducing 
the binding energies and the scattering lengths in Tables 6 and 7, we search values of the 
parameters, in Eq. (13). In doing so, we assume that the first term (Vi of Eq. (13)) is an 
attractive potential like OBEP and the second term (V2) is repulsive like the QCM repulsion. 
Accordingly, 62 is taken to be the radius of the quark wave function, 62 = 0.5fm. On the other 
hand, 61 = 0.9 fm is chosen to represent a typical range of the meson exchange potential 
which couples T,cN and Y,*N to AcN. 

The remaining parameters of the effective potential are chosen as 


= -150.0[MeV], = 109.0[MeV], V(+ = -149.0[MeV], V2+ = 98.5[MeV]. (16) 


aO-T 


1+ 


A+ 


Fig. 9 illustrates the shape of the effective potentials. 
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Fig. 10 Jacobi coordinates for AcNN system. 


Alternatively, one can express the effective potential with the strengths of the spin- 
independent terms and the spin-spin terms, 

1/^ = -149.25[MeV], = 0.25[MeV], = 101.125[MeV], = -2.625[MeV]. (17) 

It is obvious that the spin dependent potential is weak. This is a consequence of the heavy 
quark spin symmetry [20-23]. 


4-3. AcNN bound states 

Next, we study the AcNN three-body system by using the above AcN effective potential. In 
the calculation, we adopt the Jacobi coordinates shown in Fig. 10. The angular momenta, I 
and L are defined corresponding to r and R, respectively. On the other hand, we define the 
spin for two nucleons as Snn, and the total spin of the systems as Stot- For the isoscalar 
case, the binding energy is calculated from the threshold of Ac plus the deuteron, while for 
the case 1 = 1, the calculation is from the threshold of the bound AcN (binding energy=1.56 
MeV) plus a nucleon. 

In this study, we assume that the orbital angular momenta, I and L, are both zero. Then 
Snn and I are related and possible are 


1 = 0 
I = 1 


Snn = 1 , = k 


Snn = 0 , J'^ = 


2 T 


(18) 


There are various potential models for the NN system [15-19]. Here we employ the Min¬ 
nesota potential [15]. It describes the deuteron only with S-wave component, as the D-wave 
contribution is effectively included in the central potential. 



{l + P[^)Vt + 



(19) 

1 = VoRe-^^^\ Vtir) 

= Vote 


(20) 

For = 

= 200.0 [MeV], 

HR = 

1.487[fm-2], 


Vot = 

-178.0[MeV] 

, Kt = 

O 

CO 

CO 

r 

1 

to 


b/,. = 

-91.85[MeV] 

, Ks = 

: 0.465[fm"2]. 

(21) 
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r 


A^NNil = 0) 


1 r 


w/o Coulomb w/ Coulomb A^nn 
A^ + N + N 


A,NN(I = 1) -1 

A^np A^pp 




Ac +Deuteron -2.22 


> 

0) 


-5 - 


B -10 
0) 
c 
LU 

CTi 

^ -15 
m 


-20 


-21.78 V,+ 


-22.46 


-25 


-20.43 

-21.09 


AcU + N -1.56 


-10.26 V-+ 


-9.10 V, 


A^p-l-iV -0.72 


-7.17 7:+ 


Fig. 11 AcNN binding energies. The dotted lines are the Ac-deuteron threshold that is 
2.22 MeV below the AcNN threshold, the threshold of the (Acn) I"*" bound state plus a 
nucleon, i.e., 1.56 MeV below AcNN, and the threshold of the (Acp) I"*" bound state plus a 
nucleon, i.e., 0.72 MeV below AcNN 


We again use GEM [14] to solve the three-body problem. The binding energies of the 
AcNN three-body system with the AcN effective potential are shown in Fig. 11. The left 
two plots show the results for / = 0 without and with the Coulomb potential, respectively. 
The others represent the 1 = 1 results for AcUn, AcUp, and Acpp. For each system, we find 
a bound state. Because the spin-spin interaction between Ac and N is weak ( Table 7 ), 
the binding energies of = 2 and | are close to each other for the case of / = 0. These 
results agree well with the heavy quark spin symmetry [20]. In the present calculation, we 
find that the | state is lower in energy than the ^ state by about 0.7 MeV. These results 
are consistent with a recent three-body calculation [10], which also favors the I"*" AcNN as 
the lowest state. Effect of the Coulomb potential between Ac and proton is smaller than the 
difference between the 1 = 0 and 1 results. 

Tables 9 and 10 show the mean distances in the bound state, where r is the root mean 
square (rms) distance between the nucleons, and R is rms distance between Ac and the center 
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Binding energy [MeV] 

(From threshold) 

r [fm] 

R [fm] 

AcUp w/o Coulomb 

~ 2 

21.78 

19.56 

1.91 

1.34 

AcUp w/o Coulomb 

3 + 
~ 2 

22.46 

20.24 

1.90 

1.32 

AcUp w/ Coulomb 

1 + 

2 

20.43 

18.21 

1.93 

1.36 

AcUp w/ Coulomb J'" 

_ 3 + 

2 

21.09 

18.87 

1.91 

1.34 


Table 9 AcNN binding energies and particle distances for the 1 = 0 bound states 



Binding energy [MeV] 

(From threshold) 

r [fm] 

R [fm] 

AcUn J'^ 

2 

10.26 

8.70 

2.62 

1.64 

Acup 

1 + 

~ 2 

9.10 

7.54 

2.67 

1.68 

AcPP J'^ 

1 + 

2 

7.17 

6.35 

2.78 

1.75 


Table 10 AcNN binding energies and particle distances for the 1 = 1 bound states 


of NN. These Tables illustrate that the mean NN separations, r’s, are signihcantly smaller 
than the mean distance of p — n in the deuteron, about 3.8 fm. Namely, the Ac attraction 
to the nucleon makes the AcNN system shrink. It is also observed that R is smaller than r. 
Thus we can draw an intuitive picture that the nucleons go around Ac sitting at the center in 
the AcNN system. This property is also seen in the hypernuclear system, but the attractive 
force in the charmed nuclear system is stronger than that in the hypernuclear system. 

We may check whether these results depend on the choice of the NN potential. Here 
we have used the Minnesota potential, which is known to have a weak repulsive core. For 
comparison, we replace the Minnesota potential by the AV8' potential [24] and apply it to 
the 1 = 1 Acun system. For the 7 = 0 AcNN systems, the comparison is not appropriate 
because the AV8' contains the tensor force which make the D-wave components couple. 

Then we obtain the binding energy BE= 11.93 [MeV] and the rms distances as r = 2.41 
[fm] and R = 1.54 [fm]. These results can be compared to the first line of Table 10. One sees 
that they are consistent with each other. Thus we conclude that the qualitative behavior of 
the AcNN bound state is independent of the choice of the NN potential. 

5. Conclusion 

We have examined the interactions between the ground-state charmed baryons Ac, Sc and S* 
and the nucleon N. Potential models are proposed, which are composed of the long-range 
one-boson exchange (OBE) force and the short-range quark exchange force based on the 
quark cluster model (QCM). We also include the Coulomb potential between the charged 
baryons. The parameters in the model, the cutoffs and the coupling constant of a meson, are 
determined so that the NN data are reproduced in the same model. By fitting the deuteron 
binding energy and the S-wave scattering lengths, four parameter sets are obtained. 

We have applied these potentials to the two-body AcN — JlcN — S^A^ coupled system. 
The coupled-channel Schrodinger equation is solved to a good precision by the Gaussian 
expansion method. We find shallow bound AcN states both in = 0'*' and I"*" in two 
of the four parameter sets. Those parameters are the cases that the ranges of the QCM 
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repulsion are small and thus the quark exchange effect is weaker than the other two. The 
difference between the J'^ = 0^ {AcNS q)) and T*" {A^NSi)) systems is small, which is a 
consequence of the heavy quark spin symmetry. 

It is not surprising to have a bound state in these channels, because the YcN interactions 
are similar to those oi YN except for the K exchange part, while the kinetic energy is 
suppressed by the heavier mass. Furthermore, the effect of the channel coupling to T,*N 
is significant due to the strong tensor force coming from the one-pion exchange. We then 
conclude that shallow A(.N bound states may exist. 

Encouraged by the possible existence of the two-body bound states, we further consider 
three-body A^NN bound states. This is the lightest “charmed nucleus”, which corresponds 
to the hyper-triton {Apn) bound state in the strangeness sector. In order to simplify the 
calculation, in the present approach, we first construct an effective one-channel {AcN) poten¬ 
tial from the fully coupled {AcN — YcN — Y*N) calculation. Using the effective one-channel 
potential, we have solved the three-body Schrodinger equation with the GEM and have 
obtained bound states with the binding energy about 20-22 MeV for 1 = 0, and 7-10 MeV 
for 1 = 1 from the AcNN threshold. The corresponding wave functions indicate that the Ac 
baryon makes the size of the NN system significantly smaller owing to attractive force. In 
order to confirm these results, further studies with the full and explicit couplings with the 
YcN and Y^^N channels will be necessary. 

Acknowledgment 

S. M. is supported by the RIKEN Junior Research Associate Program. This work is sup¬ 
ported in part by JSPS KAKENHI Grant Nos. 25247036, and 24250294. Y. R. L. was 
supported by the “Program for Promoting the Enhancement of Research Universities ”at 
Tokyo Institute of Technology in 2014 and partly by NNSFC(No. 11275115). This work is 
also supported by the Research Abroad and Invitational Program for the Promotion of Inter¬ 
national Joint Research, Category (C) and the International Physics Leadership Program 
at Tokyo Tech. 


14/26 



References 

[1] D. Acosta et al. (CDF II Collaboration), Phys. Rev. Lett. 93, 072001 (2004). 

[2] A. V. Manohar and M. B. Wise, Nucl. Phys. B 399 17 (1993) 

[3] S. Ohkoda, Y. Yamaguchi, S. Yasui, K. Sudoh, and A. Hosaka, Phys. Rev. D 86, 034019 (2012) 

[4] C.B. Dover and S.H. Kahana, Phys. Rev. Lett. 39, 1506 (1977). 

[5] G. Bhamathi, Phys. Rev. C 24, 181 (1981). 

[6] H. Bando and M. Bando, Phys. Lett. B 109 , 164 (1982); H. Bando and S. Nagata, Prog. Theor. Phys. 
69 , 557 (1983); H. Bando, Prog. Theor. Phys. Snppl. 81 , 197 (1985). 

[7] B.F. Gibson, C.B. Dover, G. Bhamathi, and D.R. Lehman, Phys. Rev. G 27 , 2085 (1983). 

[8] Y.A. Batusov et al., JETP Lett. 33, 56 (1981). 

[9] Y. R. Lin, and M. Oka, Phys. Rev. D 85, 014015 (2012). 

[10] H. Garcilazo, A. Valcarce, and T.F. Garames, Phys. Rev. G 92 024006 (2015). 

[11] A. Gal, H. Garcilazo, A. Valcarce, and T. Fernandez-Carames, Phys. Rev. D 90, 014019, (2014) 

[12] M. Oka, and K. Yazaki, Phys. Lett. B 90, 41 (1980); M. Oka, K. Shimizn, and K. Yazaki, Nucl. Phys. 
A 464 700 (1987); M. Oka, K. Shimizu, and K. Yazaki, Prog. Theor. Phys. Snppl. 137 1 (2000); M. 
Oka, Nncl. Phys. A, 881, 6 (2012) 

[13] S. Takenchi, O. Morimatsn, Y. Tani and M. Oka, Prog. Theor. Phys. Snppl. 137, 83 (2000). 

[14] E. Hiyama, Y. Kino, and M. Kamimnra, Prog. Part. Nncl. Phys. 51 , 223 (2003). 

[15] D. R. Thompson, M. Lemere, and Y. C. Tang, Nnci. Phys. A 286, 53 (1977). 

[16] T. A. Rijken, M. M. Nagels and Y. Yamamoto, Prog. Theor. Phys. Snppl. 185, 14 (2010). 

[17] R. B. Wiringa, V. G. J. Stoks and R. Schiavilla, Phys. Rev. G 51, 38 (1995). 

[18] R. Machleidt, Phys. Rev. G 63, 024001 (2001). 

[19] R. Tamagaki, Prog. Theor. Phys. 39, 91 (1967). 

[20] N. Isgur and M.B. Wise, Phys. Lett. B 232, 113 (1989). 

[21] M. Oka, Y. R. Lin, and W. Meguro, Few Body Syst. 54, 1255 (2013). 

[22] S. Yasui, K. Sudoh, Y. Yamaguchi, S. Ohkoda, A. Hosaka, and T. Hyodo, Phys. Lett. B 727 185 (2013). 

[23] Y. Yamaguchi, S. Ohkoda, A. Hosaka, T. Hyodo, and S. Yasui, Phys. Rev. D 91 034034 (2015). 

[24] B. S. Pudliner, and V. R. Pandharipande, Phys. Rev. G 56, 1720 (1997). 

[25] K. Nakamura et al. (Particle Data Group), J. Phys. G 37 075021 (2010). 

[26] A. Manohar and H. Georgi, Nucl. Phys. B 234, 189 (1984). 

[27] D. Jido, M. Oka, and A. Hosaka, Prog. Theor. Phys. 106, 873 (2001). 

[28] S. Ishida et ah. Prog. Thepr. Phys. 95 , 745 (1996). 

[29] K. Igi and K.I. Hisaka, Phys. Rev. D 59, 034005 (1999). 

[30] G. Golangelo, J. Gasser, and H. Leutwyler, Nucl. Phys. B 603 , 125 (2001); 1. Caprini, G. Colangelo, 
and H. Leutwyler, Phys. Rev. Lett. 96 , 132001 (2006). 


15/26 



Appendix A. Spin Matrix Elements of YcN OBE potential 

In this Appendix, we present the definitions of the spin dependent operators in the potentials. 
O spin is the spin-spin operator between and N: 

Ospin = ■ <T2, (A1) 


where <T2 is the Pauli matrix of the nucleon spin, and Oi is the spin operator of the charmed 
baryon. 



0-1 

Si 

Si 


for Ac and Sc 

for the transition from Ac and Sc to S* 
for S* 


(A2) 


The transition spin S is defined as = S$, where is the Rarita Schwinger field, and 
is the spin wave functions of S*, 



/ 1 \ 


/ 0 \ 


/ 0 \ 


/ 0 \ 

d>(3/2) = 

0 

0 

, cI>(l/2) = 

1 

0 

, cI.(-l/2) = 

0 

1 

, $(-3/2) = 

0 

0 






\^) 


11 / 


(A3) 


Then we calculate the transition spin explicitly. 


St = -^(^4+ ist)+^(^st-ist)+st 

0 00\/oy| 0 0 0 o\ (A4) 

\0\/iooy\^o 0 \/ioyvoo o i ) 


— —/4x4j StS — —212x2- (A5) 

With e(A =-hi) =-^(l,i,0), e(A =-1) = ^(1,-i, 0), e(A = 0) = (0,0,1), and S\^ = 
(o,5t), one has 



/ 

1 

0 

\ 


/ 

0 

0 

\ 


/ 

0 

0 

\ 



0 

/T 




0 

0 





0 


S(l,+1) = 


V 3 


, S(1,-1) = 


r[ 

0 


, S(1,0) = 


V ^ 




0 

0 



V 3 




0 

/ ^ 
V 3 



1 

0 

0 

/ 


1 

0 

1 

/ 


1 

0 

0 

/ 


(A6) 


Next, we define the spin operator of S*, 

s = = {s\ycT{Sty, 

s{K) = 


(A7) 
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With Eq. (A 6 ), the explicit matrices are 


+1) — (^ 3 ; + ~ 


I 0 0 ^ 

0 0^0 


V 


crr^(l, -1) = ^ (Sa; - i'Sy) = - 


0 0 0 
0 0 0 0 
/ 0 0 0 0 \ 
0 0 0 
0 0 

0 0 

/ 1 0 0 0 \ 


2 

3 

D 2\/2 

^ 3 


V 


^ 0 
3 ^ 


<T,.^(1,0) = = - 


0 


0 0 


0 0 


0 


\ 0 0 0 -1 / 

The tensor operators can be defined similarly as follows 

3(<ti • r)(<T 2 • r) 


In An 


In Ac and Sc channels : Oten = 
S! and Sc 


— CTi • <T2, 


, 1 3 (S-r)(cr2-r) - 

Sc channels : C'ten = - 5 -S • <T 2 , 


In S* diagonal channels : Oten = 
The spin-orbit operator Ols is defined as 

Ols = L ■ 0-2. 


3(S • r)(<T 2 • r) 


- S • 0 - 2 , 


LOi is not included in the potential for this calculation. 


(AS) 


(A9) 

(AlO) 
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Appendix A.l. I = J'^ = 0+ coupled system 

Tables A1-A3 give the matrix elements of Ogpim ^tem and Ols for the channels in / = 1/2, 
= 0^. We label the relevant channels by i and j, and tabulate the ij component of the 
matrix elements (O)--. 

\ I IJ 


j 

i 

AcN{^So) 

SciV(i5o) 


KN{^So) 

-3 

-3 

0 

SciVC^o) 

-3 

-3 

0 


0 

0 

1 


Table A1 The matrix elements of the spin-spin operators {Ospin)ij for the / = = 0^ 

coupled system. 


j 

i 

AcNi^So) 

^cNi^So) 

^*cN{^Do) 

KN{^So) 

0 

0 

-V6 

ScA(i5o) 

0 

0 

-V6 

^*cN{^Do) 

-^/6 

-^/6 

-2 


Table A2 The matrix elements of the tensor operators {Oten)ij for the / = = 0^ 

coupled system. 


j 

i 

AcN^So) 

^cN^So) 

S*A(5z1o) 

AcNCSo) 

0 

0 

0 

^cNCSo) 

0 

0 

0 

^*cN{^Do) 

0 

0 

-3 


Table A3 The matrix elements of the orbital-spin operators {OLs)ij fo^’ ^ — 

0 "'“ coupled system. 
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Appendix A.2. / = = 1+ coupled system 


For / = 1/2, = 1+, 

the matrix elements 

are given in 

Tables A4-A6. 


j 

i 

AcNi^Si) 

ScA^(^Si) 

^*cN{^Si) 

AcNi^Di) 

EcNi^Di) S/iV(3Di) 


AcAf(3Si) 

1 

1 

Vi 

0 

0 0 

0 

T.cN{^Si) 

1 

1 

/s 

V 3 

0 

0 0 

0 

E/iV(3Si) 

fs 

V 3 

fs 

V 3 

5 

3 

0 

0 0 

0 

t^cNi^Di) 

0 

0 

0 

1 

1 

^ y 3 

0 

T.cN{^Di) 

0 

0 

0 

1 

^ -y s 

0 

^*cN{^Di) 

0 

0 

0 

fs 

V 3 

fs 5 

V 3 3 

0 

KNi^Di) 

0 

0 

0 

0 

0 0 

1 


Table A4 The matrix elements of the spin-spin operators {Ospin)ij for the / = = 1"*“ 

coupled system. 


j 

i 

AcNi^Si) 

ScAf(^Si) 

^*cN{^Si) 

AcNi^Di) 

EciV^Di) 


^*cN{^Di) 

AcN(-^SA 

0 

0 

0 

78 

78 

1 

73 

EciV(^Si) 

0 

0 

0 

78 

78 

\ 

V3 

73 


0 

0 

0 

1 

V3 

1 

a/3 

a/2 

3 

-72 

AcN{^Di) 

V8 

78 

1 

V3 

-2 

-2 

1 

a/6 

/? 

V 2 

^cNi^Di) 

V8 

78 

1 

-2 

-2 

1 

a/6 

fs 

V 2 

KNi^Di) 

1 

a/3 

1 

a/3 

a/2 

3 

1 

V6 

1 

a/6 

1 

3 

-1 

E*N{^Di) 

73 

73 

-72 

/J 

V 2 

/? 

V 2 

-1 

-1 


Table A5 The matrix elements of the tensor operators {Oten)ij for the I = J'^ = 

coupled system. 


j 

i 

AcNAsa 

t.cNAsa 

^Ni^SA 

AciV(3Di) 

EciV^Di) 

^*cN{^DA 

^*cNAda 

AcNi^SA 

0 

0 

0 

0 

0 

0 

0 

EciV^Si) 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

AcN^Di) 

0 

0 

0 

-3 

-3 

0 

0 

EciV^Di) 

0 

0 

0 

-3 

-3 

0 

0 

KNi^DA 

0 

0 

0 

0 

0 

3 

2 

3 

2 

^*cNAda 

0 

0 

0 

0 

0 

3 

2 

5 

2 


Table A6 The matrix elements of the orbital-spin operators {OLs)ij fo'^ ^ ~ ~ 

T^ coupled system. 
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Appendix B. Definitions of the radial functions and coupling constants 

Appendix B.l. The explicit forms of the radial functions 
The Yukawa potential functions in Eq. (1) are defined as 

p-X 

Y{x) = —, 

+ 

+ ^ + 

Yi(m, A,r) = Y{mr) - Y{Kr) - 

Y 3 (m, A, r) = Y{mr) - f—^ Y(Ar) - 

\m J 2m^ 

Z 3 (m, A,r) = Z{mr) - Z(Ar) - Y(Ar), 


H^{m, A,r) = H{mr) — ( — H{Ar) — 


A\^ . (A^ -m^ )A 

2m^ 


m 


Y(Ar)YA^,-^^. ,B1) 


In Eq. (Bl), A is a cutoff parameter introduced in the monopole type form factor 


FiQ) = 


A 2 - 


m 


(B2) 


A 2 _ g 2 ’ 

where m is the mass of the exchanged meson and q is the 4-dimensional momentum of 
the meson. We use the following values of the meson masses: mvr = 137.27[MeV], = 

600.0[MeV]. 


Appendix B.2. The coupling constants for the meson-exchange potentials 
The coupling constants are given in Tables Bl and B2, 



AcN 

YcN 

Y*N 

AcN 

0 

{-^92gA) 

-i^949A) 

YcN 


{-gi9A) 

{-939a) 




{g59A) 


Table Bl C-^ for the YcN systems 


a 

AcN YcN 

YIN 

AcN 

{2lBhcr) 


YcN 

i-lsha) 


Y^N 


i-lsha) 

Table B2 for the YcN systems 


with the following values [25-30]: and h^ is a free parameter to be determined in sect. 

g2 = -0.598 , 54 = 0.999 , gi = ^54 , 93 = \J\gi , 95 = , 

gA = 1-25 , Ib = -3.1 , Is = -21b, 


2.4 The explicit values of the coupling constants are given in Table B3 and B4 For the 
NN potential, we note that the function Y^{m, A,r) is used in place of Yi{m, A,r) in the vr 
exchange potential and we adopt the coupling constants given in Tables B5 and B6. 
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C. 

AcN 

YcN 

YIN 

Co 

AcN 

YcN 

Y*N 

AcN 

0 

0.92 

-1.53 

AcN 

(-6.2/io) 



YcN 


-1.18 

-1.02 

YcN 


(-6.2/io) 


YIN 



-1.77 

YIN 



(-6.2/io) 


Table B3 YcN potential Ctt values Table B4 YcN potential Co- values 


c. 

NN{I = 0) NN{I = 1) 

Co 

NN{I = 0) NN{I = 1) 

NN{I = 0) 

29 a 

NN{I = 0) 

-K 

NN{I = 1) 

29 a 

NN{I = 1) 

-hi 

Table B5 

Ctt for the NN systems 

Table B6 

Co for the NN systems 
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Appendix C. YcN potentials 

Here we show the shapes of the Y^N potentials in the present study. The potentials for 
= 0+ are given in Figs. 3 - 8 in sect. 2.4. 

Appendix C.l. Components of the Y^n CTNN-a potential (J'^ = 0"'"j 
The individual contributions of one-pion exchange, one cr exchange, and QCM repulsion 
are given in Figs. Cl - C3 for J'^ = 0'*' CTNN-a potential. One sees that one a exchange 
is strongly attractive, so that even with the QCM repulsion the total potential becomes 
attractive. The similar behavior is also observed for CTNN-b,c,d potentials. 



Fig. Cl TcA^-CTNN potential for A^N Fig. C2 y^iV-CTNN potential for 
single channel single channel 



Fig. C3 YcN-CTNlA potential for single channel 
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Appendix C.2. YcN potentials(J'^ = ) 

Figs. C4 - C29 show the potentials. Four lines correspond to the four choices 

of the parameter sets, a ~ d. The potential for the diagonal KcNi^Si) channel is not shown 
here as it coincides with that for A.c,N{^Sq). The off-diagonal potential AcNi^Si —^ Di) is 
zero. 
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V[Mey] , , , , , y[MeV] , , , y[Mey] V[MeV] 



Fig. C22 Sc*iV(35i)-S*A^(3L)i). 


Fig. C23 Sc*A^(35i)-S*A^(5Di). 



Fig. C24 AcN{^Di)-T.cN{^Di). 


Fig. C25 AcN{^Di)-J:iN{^Di). 


































